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1 Introduction 



Interaction between algebraic geometry of the compact Riemann surfaces and the the- 
ory of integrable systems represent nowadays a well established paradigma of modern 
mathematical physics. The analysis on the compact Riemann surfaces was essentially 
completed in 19th century in the famous classical works by Gauss, Euler, Riemann, Ja- 
cobi, Weierstrass etc. The penetration of the tools provided by their methods in the 
theory of integrable systems started in the same time. Legrange first solved the equa- 
tions of motion of the Euler top by means of the use of the Jacobi elliptic functions. Next 
essential breakthrough was achieved in 1889 by Kovalevski jjj who found the new inte- 
grable case of the motion of rigid body solvable in terms of genus 2 algebraic functions. 
Slightly later Dobriner published the explicit solution of sine-Gordon equation expressed 
in terms of two-dimensional theta-functions ||. Then Carl Neumann solved the equa- 
tions of geodesic motion on three-axis ellipsoid in terms of two-dimensional hyperelliptic 
theta-functions ||. 

The same time it started the development of the spectral theory of the Sturm-Liouville 
operators with periodic coefficients. First important steps in the field were done by Floke 
@, Lyapunov || and others (see related references in the book ||. 

In 1919th french matematician Jule Drach wrote the remarkable article (absolutely 
forgotten for next 60 years) devoted to the construction of explicitly solvable Sturm- 
Liouville equations associated to hyperelliptic curves of any genus @. 

One of the subsequent developments of the discovery of the inverse scattering trans- 
form (1ST) method by Gardner, Green, Kruskal and Miura (8| was the complete under- 
standing of the deep interplay between the classical works listed above. This connection 
now became a part of what is called algebrogeometric approach to the solution of nonlin- 
ear differential equations. This field of activity was initiated in 1974-1976 by the works 
of (in chronological order) Novikov, Lax, M.Kac, Its, Matveev, Dubrovin, McKean, van 
Moerbeke and Krichever. 

Speaking about the aspects of the theory connected to the explicit solutions of these 
equations one should mention the formulas for solutions of the Korteveg-deVries equa- 
tion with periodic initial data obtained in |J . Soon after the same kind of formulas were 
derived for numerous other equations integrable via 1ST method, including Non- linear 
Schrodinger [pX| ], Sine-Gordon pi]] , Toda lattice |L|, Kadomtzev-Petviashvili [13| and so 
on. The generic feature of the algebro-geometric solutions of the systems listed above 
is that they are parametrised by the fixed algebraic curves and the associated dynamics 
is linear on their Jacobians. These results found many beautiful and unexpected appli- 
cations in various branches of modern mathematics and physics, including differential 



geometry of surfaces [14] and algebraic geometry (Novikov hypothesis |15| ). 

It is also relevant to notice that simultaneously with the soliton theory, the algebro- 
geometric methods were actively applied in the framework of twistor theory (which, 
essentially, also deals with integrable equations, and exploits the methods very closely 
related to ones used by integrable systems community, but essentially focuses on the 
global properties of the solutions). Probably, the main achievments in this direction 
were classifications of the instanton and self-dual monopole configurations (see for more 
details the book and references therein). 

Technically and conceptually new development was the application of algebrogeomet- 
ric ideology to the Einstein equations of general relativity in presence of two commuting 
Killing vectors (Ernst equation) . The embedding of this equation in the framework of the 
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1ST approach was initiated by Belinskii-Zakharov Jl7| and Maison p^J Q. The character- 
istic feature of the associated zero-curvature representation is the non-trivial dependence 
of the associated connection on the spectral parameter. Namely, the connection lives on 
the genus zero algebraic curve depending on space-time variables. This peculiarity en- 
tails the drastic change of the construction and qualitative properties of algebro-geometric 



solutions first obtained in 1988 21]. Dynamics in these solutions is generated by hy- 
perelliptic curves with two coordinate-dependent branch points. The algebro-geometric 
solutions of the Ernst equation do not possess any periodicity properties which were in- 
evitable for all KdV-like cases studied before. Moreover, we can explicitly incorporate in 
the construction an arbitrary functional parameter which never appears in the traditional 
KdV-like 1+1 integrable systems. Wide subclass of the obtained solutions turns out to 



be asymptotically flat [22]. As a simple degenerate case the algebro-geometric solutions 



contain the whole class of multisoliton solutions found by Belinskii and Zakharov [17] 



Simplest elliptic solutions were studied in [22]. Despite many interesting properties, 
they contain ring-like naked singularities making it difficult to exploit them in a real 
physical context. More realistic physical application of the particular family of the class 
of algebro-geometric solutions came out from the series of papers of Meinel, Neugebauer 
and their collaborators starting from 1993 [p3| , p3f| . These works were devoted to the 
investigation of the boundary-value problem corresponding to the infinitely thin rigid 
relativistically-rotating dust disc. The explicit embedding of the dust disc solution in the 
formulas of pO] was given in pq]. In the subsequent series of papers Klein and Richter 



(see, for example, [26, |27|]) established the link of the algebrogeometric solutions of Ernst 
equation given in p0| ; p5f | with the scalar Riemann-Hilbert problem on hyperelliptic curve 
and further discussed potential applications to rotating bodies. 

Purpose of this paper is twofold. First, we extend the link between isomonodromic 



solutions of the Ernst equation and classical Schlesinger system [28|, established in the 



paper |3C], on the level of algebro-geometric solutions. Namely, we show how all algebro- 
geometric solutions of the Ernst equation may be obtained from the algebro-geometric 
solutions of the Schlesinger system found in the recent paper p9|. This allowes to get 



remarkably short expression for algebro-geometric solutions of the Ernst equation: 



*(f,o = — r-b<-' (L1) 

in terms of theta-functions with constant characteristics r, s 6 R 90 , associated to hyper- 
elliptic curve Cq of genus go with two coordinate-dependent branch points: 

2<7o 

v 2 = ( w - £)( w - $Y[( W - Wj ) . (1.2) 

i=i 

Applying certain limiting procedure to this solution, we can get general algebro-geometric 
solution of the Ernst equation. 

Second, we give explicit expressions for all metric coefficients corresponding to algebro- 
geometric solutions of the Ernst equation. The most non-trivial part is the calculation 



3 Coming back to the works of classics of differential geometry of 19th century it is curious to mention 
that essentially the same Lax pair appeared in the work of Bianchi Jl!^ in the study of so-called Bianchi 
congruences. 
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of so-called conformal factor; for this we use the link between the conformal factor and 
tau- function of the Schlesinger system, established in |3(J, and the formula for the tau- 
function obtained in [29|. The final formula for the conformal factor corresponding to 



solution fll.lD looks as follows: 



e [g] (o[B) ^ 

\J det Aq 



e 2k = ^Lviv^j yi \ Wj - r 1/4 , (i.3) 

3=1 



where Ao - is the matrix of a-periods of holomrphic differentials w 3 dw/v , j = 1, ... ,go 
on Cq, and the theta-function is associated to the hyperelliptic curve C of genus 2<?o — 1 
defined by the equation 

29+2 



y 2 = ]1 (7-7j) 



3=1 

where 



Appropriate limiting procedure allows to deduce from this formula the expression for 
conformal factor corresponding to general algebro-geometric solutions. Notice, that it 
seems to be impossible to give simple expression for the conformal factor entirely in terms 
of the objects corresponding to curve Co, which was the obstacle for deriving this formula 
without understanding the link between the conformal factor and the tau-function. 

2 Schlesinger system and stationary axisymmetric Einstein 
equations 

2.1 Schlesinger system 

Consider the following linear differential equation for function ^(A) E SL(2, C): 

f- = A( 7 )« (2.1) 

where 

A A- 

Ml) = E^ 1 - (2-2) 
f^i 7 - 7j 

and matrices Aj G sl(2, C) are independent of 7. Let us impose the initial condition 

*( 7 = oo)=J (2.3) 



Function ^(7) defined by ( p.l| ) and (^^) lives on the universal covering X of CP 1 \ 
{71, . . . ,7tv}. The asymptotical expansion of ^(7) near singularities jj is given by 

9 (7) = Qi(/ + 0(7 - 7i))(7 " Ijf'Cj (2.4) 
where Qj, Cj G SL(2, C) and Tj is traceless diagonal matrix. Matrices 

M j = Cr 1 ^C j , j = l,...,N (2.5) 
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are called monodromy matrices. 

The assumtion of independence of all matrices Mj of the parameters 7j : 

dMi , 
= 2.6 

d-fk 

is called the isomonodromy condition; it implies the following dependence of VP (7) on jj, 
which can be deduced from (2.4): 

d^/ A 



j 



djj 7 - 7j 



VP (2.7) 



The compatibility condition of (^J) and (2.7) is equivalent to the Schlesinger system [E8| 
for the residues Af. 

dAj _ [Aj,Aj] ^ 8Aj _ y [Ai,Aj] ^ 

The Schlesinger system admits the "multi-time" Hamiltonian formulation ( |p2[ ) with 
respect to the following Poisson structure: 

{A],A b k } = ffA%8 jk (2.9) 

where f® b are structure constants of sl(2). Evolution with respect to "times" jj is 
described by the Hamiltonians 

tvAjA k 



The function t({7,-}), generating Hamiltonians iifj according to equations 

^-Inr = Hj (2.11) 



is called the r- function of Schlesinger system. Compatibility of equations (2.11) follows 
from the Poisson commutativity of all Hamiltonians Hj. 

2.2 Stationary axisymmetric Einstein equations from Schlesinger sys- 
tem 

The Einstein equations for the line element 

ds 2 = f- l [e 2k {dz 2 + dp 2 ) + p 2 d<p 2 ] - f(dt + Fdip) 2 (2.12) 

where all metric coefficients /, k, F are assumed to depend only on p and z, reduce to 
the Ernst equation 

1 

—1 

P 

The metric coefficients may be restored from the complex-valued Ernst potential £(z,p) 
according to the following equations: 



{£ + £")(£„ + -E, + £„) = 2(£j + £j). (2.13) 
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where £ = z + ip. 



Equation ( 2.13 ) is the compatibility condition of the following linear system: 

-5 qi- = —S 

1-7 4 1+7 

where 



1^ = ^ * *| = ^ (2.15) 



7 = { A - + ^(A-OCA-I)} , (2.16) 
A G C is a spectral parameter and 

G = ITS ( i(£ - £) ^ P ) (2 ' 17) 
In terms of the matrix G the Ernst equation may be equivalently rewritten as follows: 

(pGpG-X + {pGzG-% = (2.18) 

Relationship between solutions of Schlesinger system and Ernst equation was revealed 

in 



Theorem 2.1 Let {Aj} be some solution of the Schlesinger system (2.8) and ^(7) be 



related solution of equation (2.1) satisfying the following conditions: 



**(-)^(0)" 1 ^(7) = / , (2.19) 

7 



*(-7) = *(7) • (2.20) 

Let in addition 7^ = £, £), Xj G C /or aZZ j, w'i/t function 7(A, £, £) given by ( \2.1&\ ). 
Then 

G(e,o^*(7 = o,e,o (2.21) 

is a solution of Ernst equation (\2.18\ ). Function ^(j, {7j}) solves the associated linear 
system fi2.1!\) . Metric coefficient e 2k of the line element Ij2.12\ ) is related to the r -function 
of the Schlesinger system as follows: 



e 2k 



an ^ r, (2,2) 



where C is a constant of integration. 



3 Solutions of the Schlesinger system in terms of theta- 
functions 

Let TV = 27 + 2. Define hyperelliptic curve C of genus 7 by the equation 

2g+2 

11(7-7,) (3.1) 



w 2 



with basic cycles (a,-, &,-) chosen according to figure 1. 
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The basic holomorphic 1-forms on C are given by 

t^LL, k=l,...,g. (3.2) 

w 

Let us define g x g matrices of a- and 6-periods of these 1-forms by 

A kj = I B kj = I ^L. (3.3) 

J at w Jb 3 W 

Then the holomorphic 1-forms 

u 

ir 



1 9 

dU k = -Y,(A- l ) k rf- 1 d 1 (3.4) 



satisfy the normalization conditions j> a df/^ = (5^^. 

The matrices A and i3 define the symmetric g x g matrix of 6-periods of the curve C: 

B = A' X B . 

Let us now introduce the theta function with characteristic [Jj] (p £ C 9 , q € C 9 ) by the 
following series, 

6[P](z|B)= exp{7ri(B(m + p),m + p) + 2vri(z + q, m + p)}, (3.5) 

for any z E C 9 . It possesses the following periodicity properties: 

e[P](z + e,) = e 2 ^e[P](z), (3.6) 

g] (z + Be,) = e - 27r ^e- 7riB «- 27r ^e[P](z), (3.7) 



e, = (0,...,l,...,0)* (3.8) 



where 

(1 stands in the jth place). 

The theta-function with characteristics is related as follows to the theta-function 
without characteristics: 

0[P]( Z ) = 0(Z + Bp + q ) e ™<Bp,P>+2™<P,z+q> ( 39 ) 

Cut curve C along all basic cycles to get the fundamental polygon C. For any mero- 
morphic 1-form dW on C we shall use the notation 

rP 

W\n= / dW 



where the integration contour lies inside of £ (if dW is meromorphic, the value of this 
integral might also depend on the choice of integration contour inside of C). By U\q we 
shall denote vector with components Uj\n. The vector of Riemann constants correspond- 



ing to our choice of the initial point of this map reads as follows [31]: 



K = -B(ei + ... + e g ) + -(e a + 2e 2 . . . + ge g ). (3.10) 
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The characteristic with components p G C 9 /2C 9 , q G C 9 /2C 9 is called half-integer 
characteristic: the half-integer characteristics are in one-to-one correspondence with the 
half-periods Bp + q. If the scalar product 4(p, q) is odd, then the related theta function 
is odd with respect to its argument z and the characteristic [q] is called odd, and if this 
scalar product is even, then the theta function 0[q](z) is even with respect to z and the 
characteristic [q] is called even. 

The odd characteristics which will be of importance for us in the sequel correspond 
to any given subset S = {7^, . . . , 7i 9 _i} of g — 1 arbitrary non-coinciding branch points. 
The odd half-period associated to the subset S is given by 



g-i 

Bp s +q s = J2u\ 

3=1 



K. 



71 



(3.11) 



,dU g y. Denote by f2 7 C 



where dU = (dU\, . . . 
7 = 00, such that 

Let the 2x2 matrix- valued function $(7) be defined in the domain f2 7 of the first sheet 
of £ by the following formula, 



the neighbourhood of the infinite point 
7 does not overlap with projections of all basic cycles on 7-plane. 



$( 7 e 



^(7) 



tf>(-f) 



(3.12) 



where functions ip and tp are defined in the fundamental polygon C by the formulas: 



v(7) = ©KI (u\i + u\%\b)q[s\ (u\ 
Hi) = e IJ] (u\l + u\% b) e [s\ (u\ 



17 
1 7^ 



B 



B 



(3.13) 



(3.14) 



with two arbitrary (possibly {7j}-dependent) points 7^, 7^, G £ and arbitrary constant 
characteristic [q] ; * is the involution on L interchanging the sheets; 



8 [5] 



e 



(») 



where odd theta characteristic ^ s corresponds to an arbitrary subset S of g — 1 branch 
points via Eq. ( |3.11| ). 

Since domain fi 7 does not overlap with projections of all basic cycles of C on 7-plane, 
domain C* does not overlap with the boundary of £, and functions ^(7*) and ^(7*) in 
fTT2p are uniquely defined by ( gig ), ( gig ) for 7 G 7 . 

Now choose some sheet of the universal covering X, define new function ^(7) in 
subset Qry of this sheet by the formula 



r, \ /det$(oo 1 )^_ 1 . . 

^7G^ 7 ) = ^^A_Z<1> '(oo>(7) 

and extend on the rest of X by analytical continuation. 



(3.15) 



Function ^(7) ( 3.15 ) transforms as follows with respect to the tracing around basic 
cycles of L (by T a . and T&. we denote corresponding transport operators): 



r a .[*( 7 )] = *( 7 )M a • 



r 6 .[*( 7 )] = *( 7 )M 6 . 



j = 1, ... ,2g - 1 , 
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where 

^ = ( e~™* ' M » = [ A 1 (:U,i) 



The following statement was proved in the paper [29|: 

Theorem 3.1 Let p,q £ C 9 be an arbitrary set of 2g constants such that [q] is not 
half-integer characteristic. Then: 



1. Function *S?(Q G X) defined by fjS.laj) is independent ofj^ and 7^ and satisfies the 
linear system (2.1) with 

Aj = res| 7=7j {* 7 * -1 } . (3.17) 
which in turn solve the Schlesinger system 



2. Monodromies (2.5) of^{^) around points jj are given by 

— m 



where constants m,j may be expressed in terms of p and q (see j^W). 



3. The T-function, corresponding to solution (3.17) of the Schlesinger system, has the 
following form: 

r({ 7j }) = [detA]^ " ^)~ le S (°I B ) ( 3 - 19 ) 

j<k 

4 Solutions of the Ernst equation in terms of theta-functions. 
Formulas for the metric coefficients 

According to the relationship between Schlesinger system and Ernst equation given by 
the theorem |2.1| , we can derive solutions of Ernst equation in terms of theta-functions 
from construction of the theorem |3.l| . The necessary additional work to do is to choose 
the parameters of the construction (i.e. the constants Aj and vectors p, q) to provide 
the constraints (|2.ig| ) and ( |2.20D . To get these constraints fulfilled we have to assume 
that the curve £ is invariant under the holomorphic involution a acting on every sheet 
of C as 

o : 7--, (4-1) 

7 

and anti-holomorphic involution \x acting on every sheet of C as 

H : 7 -» -7 . (4.2) 

Constraints ( ^19|) and (gjg ) turn out to be compatible with each other only if the genus 
g is odd: 

g = 2g - 1 . (4.3) 
We shall enumerate the branch points 7,-, j = 1, . . . , 4<7o m such an order that 

7; = 77+2 S o> 3 = !,•••, 250 (4.4) 

and for some k < go 

7jG«, l<j<2k; 72^+1 + 72^+2 = 0, 2/c + 1 < j < 2g - 1. 
We shall now distinquish two cases: 
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1. Curve C is non-separable ^ with respect to the action of anti- involution /x, i.e. k > 1. 
Then the basic cycles (aj, bj) , j = 1, . . . , 2go — 1 on C may be chosen as shown in 
figure 2a. 

2. Curve C is separable with respect to the action of anti-involution /i, i.e. k = and, 
therefore, none of the points jj lie on the imaginary axis. In this case we choose 
the basic cycles on C as shown in figure 2b. 

In both cases the basic cycles transform in the following way under the action of the 
involution a: 

o-(ai) = -oi a{bi) = -h (4.5) 
a( aj ) = a j+go -i a(bj) = b j+go -i , 2 < j < g (4.6) 
Constraint fl2.19|) leads to the following equations for M a . and M b . ( 0.16 ): 

K M °i+«,-> = 1 ; H^Vso- = 1 . 2 < i < 50 . 

(Equations M* x = M ai and M£ = , which arise from the calculation of monodromies 
of ( ^.19|) along the basic cycles a\ and b\, are identically fulfilled) 

In turn, we get for components of the vectors p and q the following equations 

Pj + Pj+go-i = ; Qj + Qj+go-l = > 2 < 3 < 9o (4.7) 



It remains to derive constraints on p and q imposed by the "reality conditions" (|2.20|) , 
We shall consider two cases separately. 

1. Non-separable case (k > 1). The basic cycles of C shown in figure 2a behave as 
follows with respect to the action of the anti-involution ji: 

fj,(aj) = -aj , Mj ; n(b\) = h + 2a 1 

fi(bj) = bj, 2<j<k; (i(bj) = bj - ctj , k + l<j<g 

(since we assumed already the invariance of C under the involution er, it is sufficient 
to determine the action of [/, only on the first go cycles bj). Now from ( 2.20 ) we get 
the following equations for monodromies: 

M a . =M- 1 , Vj; M bl =M bl Ml; 

M bj =M bj , 2<j<k; M b . =M 6 .M-. 1 , k + 1 < j < g , 
which is equivalent to the following conditions imposed on p and q: 

Pj € R , Vi; = pi ; 

%=0, 2<j<k; Uqj = -- Pj , k + l<j<g . (4.8) 



4 Curve admitting anti-holomorphic involution is called separable, if it is divided into two pieces by 
the ovals invariant with respect to the anti-involution, and non-separable otherwise. 
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2. Separable case (k = 0). In this case we shall choose the basic cycles according to 
figure 2b. They transforms under the action of ji in the following way: 



9o 



Vj; 



K b i) = h + J2( a i+9o-i ~ 



ai 



1=2 



90 



^(bj) = bj - aj + ^ai , 2 < j < 



(Jo 



i=i 

which leads to the following conditions on p and q: 

90 



Pj e 



Vj; 



3fy 



-. go 

It, 



vi ; 



IE 



Pi 



Pj 



(4.9) 



1=2 



1=1 



Now, taking into account the theorems JO and 3.1, we get the following 



Theorem 4.1 Let the genus of curve C be odd: g = 2go — l; and the basic cycles be chosen 
according to figure 2a if C is of non-separable and figure 2b if C is of separable type. Let 
p, q € C 290 ^ 1 be arbitrary constant vectors satisfying conditions j\4- r \) . Let in addition 
p, q satisfy conditions ( \j-S^ ) if C is non-separable and if £ is separable. Define 



function ^(7) by the expressions A8.W , H3.1$ ), j[sTji) and ( \3.l3j) . Then the function 

G(x,p) = *(x,p,7 = 0) (4.10) 

satisfies the Ernst equation ( 2.1$) and may be represented in the form (2.17). 

The expression for the Ernst potential £ may be obtained from ( 4.10[) in terms of 
theta-functions associated to curve C. However, it may be essentially simplified if we 
make use of invariance of curve C under involution a and use the spectral parameter 
A from ( 2.16| ) instead of 7. Namely, curve £ may be represented as fourfold covering 
of A-plane; its Hurvitz diagram is shown in figure 3. In this realization involution a 
interchanges the sheets 1 3 and 2 «-> 4; involution * interchanges the sheets 1 «-> 2 and 
3^4. 

Now introduce new hyperelliptic curve Cq of genus go defined by equation 



230 



(A-e)(A-oII( A -^) 

3=1 



(4.11) 



Curve Cq together with canonical basis of cycles (a°,6°) is shown in figure 4. Introduce 



the dual basis dV = {dV\, . . . , dV go ) of holomorphic 1-forms on Cq by 



^ 90 

dVj = - Yji^jkX^dX , j = 1, ... )5o 



k=l 



where 



(a: 



kj 



X k ~ l dX 



j,k = 1, 



,9o , 



(4.12) 



(4.13) 



and corresponding matrix of 6-periods Bq. 
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Curve C is twofold non-ramified covering of Co'- II : C — > Cq, such that the points of 
C related by the involution a project onto the same point of Co, namely, the A-sheets 1 
and 3 of £ are projected onto the 1st sheet of Co, and sheets 2 and 4 of £ are projected 
onto the 2nd sheet of Co, preserving projections of corresponding points on A-plane. 
Anti-involution fj, inherited from C acts on every sheet of Co as A —* A. 

Existence of reduction ( 2.19| ) allows to alternatively express function ^ in terms of 
theta-functions associated to the curve Cq. Denote by ft\ the neighbourhood of the point 
A = oo being the projection of the domain fi 7 into A-plane. Define function <I?o(A G Q\) 
in the domain £1\ lying on the first sheet of Co by the following formula: 



(po(X) M**) 
V>o(A) ^ (A*) 



<MA E n A ) = 

where involution * inherited on Co from C interchanges the A-sheets of Cq', 



?o(A) = eE] (v\*\bo) 



V>o(A) = G 



VI 



¥>o(A* 



-v| 



B 



V>o(A* 



*e E] 



-v\Uv 



B 



(4.14) 



(4.15) 



(4.16) 



for A £ Slj, and constant vectors r, s G C 90 satisfy the following reality conditions: 

r G M 90 (4.17) 



% = E^> ^i^; E 



2 ' 



k + 1 < j < go . 



(4.18) 



1=1 i=i, i^j 

The basic cycles of curve Co behave as follows under the action of anti-involution \i: 

90 



-a 



o 



Vj; 



i=i 



which implies the following relations for the matrix of 6-periods of Co: 

1 



3?(B C 



1 5 H ■ , 
- + — , j > k 

2 2 ' J 



o)ji = ~2 ' J - k ; ^(Bo)iZ : 
Now the relations ( fCTTD , ( gig ) may be equivalently represented as 

5R(B r + s) = ; 

this, in particular, implies the following relation between functions ipo and ^o : 



^o(^p) = MP) 



Now define function 



(4.19) 



(4.20) 



(4.21) 



and extend it to the universal covering X by analytical continuation to get function 
tyo(P G X). The following statement takes place: 
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Theorem 4.2 Let function Vl/(P G X) 6e defined by the formulas (\3.12Q , (\3.1S\) , ( \3.14\) , 
( 3.15 ), and function ^q(P G X) be defined by formulas ( 4-14 ), ( U-15 ), ( (.21 ). Let the 
components of vectors p,q£ C 290_1 and r,s G C 90 be related as follows: 



Pi 



go 
1=1 



<ll 



-2si ; 



Pj 



si 



2<j<g (4.22) 



(relations ( \4- % ) automatically give other g — 1 components of p and q). Require in 
addition that [g] is not a half-integer characteristic. 
Then functions ^ and $o coincide: 



qj(p eX) = V (P g X) . 



(4.23) 



Proof. The functions (po and ipo were chosen in such a way that the analytical con- 
tinuation of from £l\ on the whole universal covering X does not violate reduction 
restriction ( 2.1 9| ) . Taking into account the non-triviality assumption of non-coincidence 
of characteristic [£] with any half-integer characteristic, and coincidence of the normal- 
ization conditions of and $q at 7 = 00, it is enough to show that monodromies of 
and ^0 along the first go pairs of basic cycles of curve C coincide. For our choice of the 
basic cycles on C and Co we get the following relations between the basic cycles of Co 
and C: 



U( ai ) = -{ai + ... + a u go ) , n(6 1 ) = -26^, IL{ aj ) = a" , 1%) 
which imply the following expressions for monodromies of ^0 around (aj,bj): 



h° 



M 



ai 



exp{— 2iri fja^} , = exp{47risi<73} 



M° = exp{27rr jC T3} 







M 



exp{27r(si - Sj)a 3 } , 2 < j < go 



(4.24) 



(4.25) 



coinciding with monodromies of ^ ( 3.16 ) provided conditions ( [4.22p are fulfilled (reduc- 
tion ( |2.19| ) provides coicidence of remaining monodromies around (aj,bj), go + 1 < j < 
2go). [Notice that functions (po and ipo transform in a different way (their monodromies 
differ by sign) along the basic cycles of Co, but their pullbacks on C tranform in the same 
way along the basic cycles of C] 



Now we are in position to formulate the following statement: 
Theorem 4.3 Let r, s G C 90 be arbitrary constant vectors satisfying reality conditions 



( (.11 ), ( 4-1& ). Then the following function: 



BE] 


(y\f 


B ) 




(v\f 


B ) 



(4.26) 



solves the Ernst equation ( 2. IS ). Function ^ defined equivalently by ( 3.15 ) or ( 4-21 ) 
satisfies the linear system ( 2.11 ) with matrix G(£,£) given by l \2.1ty , and 



^( 7 = 00) = G 



(4.27) 
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Proof. The non-trivial part is to check (f4.26| ). One can represent matrix \I/(A = oo 1 ), 
given by ( 4.21| ), in the form ( 2.1 7| ) with 



£ 



Vo(oo 2 ) 

which leads to (4.26) after substitution of ( [4.15|) . 

Now it arises the non-trivial problem of integration of equations ( p. 14 ) for the metric 
coefficients F and k. The following statement shows how to find the metric coefficient 
-F(£,£) corresponding to the Ernst potential ( [4.26 ): 

Theorem 4.4 The metric coefficient F(£,£) (\2.14 ) corresponding to the Ernst potential 





({■2b) is given by 



go 



3=1 



V\f\B )} 



(4.28) 



up to an arbitrary additive constant, where Ao is the matrix of a-periods of holomorphic 
1-forms ( \j.lS ); §^ denotes derivative of theta-function with respect to its jth argument. 



Proof. Consider the following simple identity 



(4.29) 



Evaluation of (12) matrix element of the right hand side of ( 4.2S| ) at 7 = oo using 
the linear system ( p. 15 ), equation for F ( |2,14j ) and normalization condition (2.3) gives 
nothing but \F%. Therefore, we can integrate equations for F to get 

F = 2($- 1 1' A -i) i2 (7 = 00) (4.30) 

Substitution of (|415|) , (|1|) and (jl2l|) into (|43(]) leads to the following expression: 

?{ — r-(A = oo)) 



where the 1-form dWo(P E Cq) is given by 



so 



dW (P) = ^2dV J (P)^-lne[l] (v 

3=1 3 

Using formula ( 4.12| ) for basic differentials dVj , we come to ( [4.28 ). 

The next theorem shows how to integrate equation (2.14) for remaining metric coef- 
ficient k: 

Theorem 4.5 The metric coefficient e 2k ^'^ , solving equation ( 2.14 ) with Ernst potential 
defined by ( \j-26 ), is given by the following expression: 



B, 



'2k 



g g] (°i B ) n ia fi-i 



/■i 



(4.31) 
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up to multiplication with an arbitrary constant, where B is the matrix of b-periods of 
curve C ( j!?. matrix Ao of a-periods of holomorphic differentials on the curve Co is 
defined by ( 4-l3j ); the constant vectors p,q£ C 290 " 1 are expressed via ( 4.2i ) in terms of 



vectors r, s G C 90 . 

Proof. One should exploit the coincidence of functions an d ^, and substitute the 
formula for r-function of the Schlesinger system ( |3.19| ) into the relation ( p. 22 ) between 



the r-function and the coefficient e 2k (in the present case tr^4 2 = 1/8). In addition 
we have to use the following relation between determinants of matrices of a-periods of 
holomorphic 1-forms of curves C and Co- 

det A = const p 9 o det A (4.32) 



where the constant is independent of (£,£)■ Proof of ( 4.32| ) may be obtained by elemen 



tary manipulations using decomposition of holomorphic 1-forms on C into combination of 
holomorphic 1-forms on Co and holomorphic 1-forms on (£, ^-independent hyperelliptic 
curve of genus go — 1 defined by the equation 

290 

5 2 = n(A-A,). (4.33) 

3=1 



To complete the proof of ( 4.31| ) it remains to make use of the folowing relations: 



290 

= {(A - 0(A - 0r 1/2 ; IIC-w - ^ = const p~ H tli^i - ^ - ■ 

j<k j=l 



Remark 4.1 Using the standard expression for Vandermonde determinant, we get the 
following formula for det .4.0 : 



det A 



o 



ai 



dpi . . . dp. 



(4.34) 



Remark 4.2 The theta-function G g] (0|B) from ( ^3l| ) may be decomposed into com- 
bination of theta-functions associated to double matrices of 6-periods 2Bo and 2Bi of 
curves (pT|) and ( p3|) (see ||). 



Remark 4.3 Half-integer r and s. 



Remark 4.4 Degeneration \2j+i, ^2j+i Kj G R of the spectral curve Co gives after 
appropriate choice of [r, s] the family of multi-Kerr-NUT solutions. In particular, to get 
Kerr-NUT soultion itself one has to take go = 2 and 

rj = -± ; Sj = -f + idj, rij = ±1 , j = 1,2 
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Then (0~26h turns into lEOfl: 



5 _ 1-T r -i _ i(ai -02) x 1 + Qiog Y 

1 + 1 - aia 2 + i(a\ + a 2 ) 1 - aia 2 -H(ai + a 2 ) 

where a,- = nje~ 2ma J ', j = 1, 2, and 
are prolate ellipsoidal coordinates. 

Remark 4.5 Elliptic case: r = - "toron" solution; r 7^ O-non-asymptotically flat. 

5 Known form of algebro-geometric solutions of Ernst equa- 
tion 

The original algebro-geometric solution of Ernst equation found in pO] looks as follows 



in notations introduced in |22], |25| : 



S 



e 


(v\f 


+ B w 


B ) 


e 


(y\f 


+ Bw 


B ) 



M r 1-1 

4exp{^|~ 2 } (5.1) 



where dW(P) is an arbitrary 1-form on Co satisfying the following conditions: 

1. dW is an arbitrary finite, infinite or continuous linear combination, with (£,£)- 
independent coefficients, of normalized (all a-periods vanish) Abelian differentials 
on Cq of the 2nd and 3rd kind with (£, ^-independent poles and singular parts. 

2. Reality condition: 



dW{P) = dW(P) , 2mB w £ R 90 . (5.2) 

Vector 2niBw is the vector of 6-periods of differential dW: 

2m(B w ) j = J> dW(P) (5.3) 
Jbj 

5.1 Role of 1-form dW 
5.1.1 Static solutions (go = 0) 

To clarify the role of differential dW we shall consider first the simplest case g = 0, when 



(5.1) becomes real: 

]n£ = W\$ (5.4) 
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which is real as a corollary of ( |5.2| ) , and the Ernst equation linearises to the Euler-Darboux 
equation: 



1 



(di + -d p + di)]n£ = 



P 



(5.5) 



The natural question which was asked in [22] is whether it is possible to represent general 
solution of (|5.5|) invariant with respect to the change of sign of p in the form (|5.4[), The 



positive answer was given in [25|. 

Namely, if we fix some domain D in (x, p)-plane, symmetric with respect to involu- 
tion p — * —p, (which may contain the infinite point with vanishing boundary condition 
ln£ — > 0), the general solution of ( |5.5[ ) in this domain may be represented by the 

x,p— >oo 

contour integral over boundary dD (which is nothing but axisymmetric version of famil- 
iar representation of arbitrary static electric field inside of a shell as an induced field of 
certain charge distribution on the shell): 



In £ 



f{n)dn 



i8d{(k-0(k-0} 1/2 ' 

with arbitrary function f(n) defined on dD and satisfying reality condition 



(5.6) 



m = /(«) . 

Solution (|5.6| ) may be represented in the form ( |5.4|) , where dW is some locally holo- 



morphic 1-form on the rational curve Co (4.11) with go = 0. For that one should take 



dW(X) 



f(K)(dW K (\)-dW K *(\))dK, 



(5.7) 



where dW K (X) is differential of the 2nd kind on Co with unique simple pole at A = n and 
the following local expansion at A = k: 



dW K (\) 



1 



+ O(l) I dX . 



A — >k \(X — K) 2 

In the present go = case we can write dW K (X) explicitly as follows: 

Y(«)Y(A)dA 



dW K (X) 



and, therefore, 



dW« = a„ (In 7 ) | 



( 7 (A) - 7 (k))5 



7=0 
7=00 



(5- 



(5.9) 



{(K _ 0(K _ e - )} i/2 • 



(5.10) 



which shows coincidence of fl5,4|) with (|5.6[). 



The partial integration in (5.7) leads to the following alternative representation of 



dW: 



where 



dW{X)=i> h(n)dW KK *(X)dn 

'dD 



h(K) 



/'(«) 



(5.11) 
(5.12) 
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and dW KK *(X) is the standard differential of the 3rd kind given by 

dW K K * (A) = ( ^ (A) - —^—) dX . (5.13) 
V7(A)-7(k) 7(A) -7 ^/c)/ 



Integral representation ( 5,11| ) of <iVF is the counterpart of the partially integrated version 
of Q: 



ln£ = 2 / /i(k) m7(ft, £, £)(fft 



Performing in ( |5.7| ) the partial integration in a different way we can represent dW as 
a contour integral over k of meromorphic 1-forms having poles of arbitrary fixed order 
at A = k, and A = k*. We see that from the local point of view the set of 1-forms dW 
allowed in fl5.4| ) is over-complete, and in order to get general local static solution it is 
sufficient to restrict ourselves by forms dW represented by ( 5.11| ). 



5.1.2 Forms dW for higher genus 

Denote by dWQn(P) the differential of the 3rd kind on Co with poles of the first order 
at Q and R and residues +1 and —1 respectively, satisfying normalization conditions 

<f dW QR {P) = , Vj . (5.14) 
J aj 

The vector of 6-periods of (IWqr: 

27ri(BQ R )j = <b dW QR 
may be expressed in the following way in terms of basic holomorphic 1-forms dVji 

(Bqr),= f dVj. (5.15) 



R 



For go > 1 it takes place the situation similar to go = 0: locally the same 1-form 
dW{P £ Co) can be represented in many different ways as contour integral of elementary 
differentials. The representation which it will be convenient to use in this paper is a 
go > 1 counterpart of ( |5.11 ): 



dW = 6 h(K)dW KK * dK . (5.16) 

JdD 

For vector of 6-periods we have 

{B w )j = 2 I h(K)Vj\UK . (5.17) 

JdD 

Now one can easily see the link to papers |27||, where it was noticed that certain combina- 



tion of components of the function ^ from (|2.1| ) corresponding to solution (5JJ, ( 5.1 



solves scalar Riemann-Hilbert problem on curve Co with contour 3D. Partial integration 



in ( 5.16 ), ( 5.17 ) gives the following relation between our function h{n) and conjugation 



function of RH problem G{k) from |27[ ]: 

d K In G(k) = -Amh{n) . (5.18) 
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The proof of ( |5.18| ) can be obtained expressing corresponding objects in terms of the 
prime-form on Co- 

Now we are going to discuss the relationship between solutions ( f4.26| ), which we de- 
rived exploiting the link between Ernst and Schlesinger equations, and previuosly known 
class of solutions of Ernst equation 

rameters; however, it turns out that in some sense both classes of solutions coincide. 



5.1). Apparently, solutions (}4.26|) contain less pa- 



5.2 Solutions (|726|) from solutions ( |5TTD 

To describe this link we shall briefly discuss how abelian differentials of the 1st kind on 
Co arise as limits of abelian differentials of the 3rd kind. 



5.2.1 Holomorphic 1-forms as limits of meromorhic 1-forms 

Cut C along basic cycles started at the same point and consider <IWqji(P) in the funda- 
mental polygon Co whose boundary consists of the cycles a° + , a° , b® + . 

Obviously, in the limit Q — > R differential dWQR should turn into certain holomorphic 
1-form. The most naive way to perform this limit is to take Q — » R inside of Co; then 
dWQR — > and Bqr — ► 0. However, this limit may also be taken in two different 
non-trivial ways: 

1. Let Q tends to some boundary point of Co, say, to a point belonging to the contour 
and let R tend simultaneously to corresponding point of the contour . Then 



the limiting procedure does not influence the normalization conditions (5.14), and 



since in the limit <1Wqr turns into holomorphic 1-form, this 1-form must vanish: 

dW QR 0. (5.19) 

However, the vector of 6-periods of dWQR does not vanish even in the limit: from 
( p. 15 ) we see that 

(B QR ) k -> -S jk . (5.20) 



2. Alternatively, we can treat the limit assuming that Q tends to some point of Co 
on the contour + , and R tends to the same point of Co belonging to the contour 
aS . Since the poles of 6Wqr meet exactly on the cycle aP, we see that the integral 
§ a dWQR does not vanish in the limit any more. Since f dWQR still vanishes for 
all k 7^ j, we conclude that in the limit dWQR becomes proportional to dVj. To 
find the coefficient of proportionality, we calculate vector of 6-periods of dWQR in 
the limit according to ( 5.15| ) and find that 



B QR -» B 0i , (5.21) 

where Boj stands for jth column of the matrix of 6-periods. Therefore, for dWQR 
itself we have 

dW QR 2ni dV j . (5.22) 
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Now choose differential dW in solution (5.1) in the following form: 

90 

dW = J2(r j dW Q . Rj -s j dWQ. R .) 

3=1 



(5.23) 



with some constant vectors r, s £ C 90 , and take the limits Qj — ► Rj Gcft, Qj — > RjGjft 
as explained in items 1 and 2 above. In this limit, according to ( 5.19 ), (|5.2C| ), ( 5.21 ) 



(|5.22|) , we get 



— > B r + s , 

90 

W\°°l -> yr,-K-l°°2 . 

loo^ / / 3 3 loo^ 



(5.24) 
(5.25) 



According to reality conditions ( |5.2j ) we should impose the following restriction on con- 
stants r and s: 

B r + sGil 



coinciding with ( [4.19 ) . Therefore, solution ( |5.1| ) turns in this limit into 

,0 ) ix- T/i-n 






[V\f +B Q r + s 


B ) 


e 


[V\f +B r + s 


B ) 



i=i 



which coincides with ( 4.26| ) if we take into account that the theta- function with char- 
acteristics in nothing but the ordinary theta-function with shifted argument multiplied 
with certain exponential factor (|3.9| ). 

Next we shall consider more non-trivial procedure of coming from ( p~26| ) to flOP, 



(5JJ). 



6 General algebro-geometric solutions of Ernst equation as 
limits of Schlesinger-related ones 



Here we shall describe the inverse procedure: how to get the class of solutions (pM|), 
(5.16) starting from solutions fl4.26j ). 



6.1 Partial degeneration of spectral curve 

Let us consider solution (|4.26|) with curve Cq substituted by the curve L\ of genus g\ 
go + n defined by the equation 



(A-O(A-e) II (^ ~~ A?) ' 

3=1 



and choose the vectors r, s G C 9o+n in the following way: 
s = ; Tj = , 1 < j < g ; r j+go 



hj e 



1 < j < n ■ 



(6.1) 



(6.2) 
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Without loss of generality we shall assume hj S [0,1]. Denote by Bi the matrix of b- 
periods of C\ and by dV\, ■ ■ ■ , dV go + n the basis of normalized holomorphic 1-forms on C\. 
Now consider the solution ( |4.26| ) constructed from these data: 



£ 



©El 


[v\f 


B1 ) 


eg] 


[v\f 


B1 ) 



and take the limit 



Kj G 



^2g +2j+l , Mg +2j+2 

Then curve C\ turns into Cq with double points at Kj ,j 
holomorphic 1-forms of C\ turns into 



(6.3) 



(6.4) 

, n. The basis of 



dVi,... ,dV„ , 



1 



2ni 



dW, 



i ■ j k\ i ■ ■ ■ i n ■ dW^n* , 



where dV\ , . . . , dV 90 is the basis of normalized holomorphic 1-forms on /2o, and dW K . K * 
are normalized 1-forms of the 3rd kind on C$ with simple poles at Kj and k* and residues 



+1 and —1 respectively. Therefore, the matrix of 6-periods of L\ in the limit fl6,4| ) may 
be described in terms of the objects associated to the curve £,q as follows: 



(Bi)jJfc = ( B o)jfc + o(l) , 1 < J, < 5f ; 
(Bi)j- fc+90 = 2^|^ + o(l) , 1 < j < g Q , 1 < k < n ; 



(Bi) 



3+so j+go 



Til 



In |A 2ff0 +2i+i - A 2so +2j+2| + 0(1) , 1 < j < n 



(6.5) 
(6.6) 

(6.7) 



( B i)fc+9o j+go = Oi 1 ) > 1 < i ^ * < n . 
Substituting these relations in the definition of theta-function we can express the value 
of Ernst potential fl6.3|) in the limit ( p^) in terms of the objects associated to the curve 



e (v\f 


+ 2£?= 




B ) 


e (v\f 


+ 2£-= 




B ) 



^exp{f>^, re *|^} 
Jo j j=l 



(6.8) 



The formula ( |4.28| ) for the metric coefficient F transforms into 



3=1 



.</<>./ 



lnG^lf + 2^/^f 

3=1 



n dW * 



3=1 



For coefficient e 2k we get from (4.31 
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e 2k = {detA)H f[l A i -^"^©(E^^ll^ " ei 



B 



xexp{27ri^^ 2 /3 J +vriB 11 (^/i i ) + ^fcW^*!^}; 



(6.10) 

3=1 3=1 3#fe, 3,fe=l 

where as before B is 2go — 1 x 2<?o — 1 matrix of 6-periods of curve C; (3j is the second 



term of asymptotical expansion of (Bi 

1 



'9+3 g+3 



as 



\+2j+i —> \+2j+2 (j — 1, • • • , n): 



( B l)a+j 9+3 - — m \^g+2j+l ~ A g+2 j+2| + &j + o(l) . 



7T/ 
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Remark 6.1 The assumption /t £ 1 , hj E M. was only made to shorten the presenta- 
tion; one can also allow in fl6.8|) the presence of conjugated pairs Kj = «j, hj = hi which 
come out if we glue together two "vertical" branch cuts. 



6.2 Continuous limit: condensation of double points 



Now we can take a continuous limit in fl6.8| ) distributing points kj over an arbitrary 
contour dD with an arbitrary (say continuous) measure h(n) satisfying reality condition 



h(n) = }i(k) 

Then (|6.8[) turns into 

e(v\f + 2$ dD h(K)v\ldK 



®{V\T +2§ QD h{K)V\1dK 



B, 



Br 



exp 



dD 



h{K)W K 



(6.11) 



coinciding with (PI) , (fTl6l) , (|5.17|) . 

Taking continuous limit in the formulas (^^) , ( |6.10| ) , we come to the following expres- 
sions for the metric coefficients e 2k and F corresponding to the Ernst potential ( |6.11 ). 



Theorem 6.1 Coefficient F of the metric ( 2.1^ ), corresponding to solution (6.1l) of th 
Ernst equation, is given by 

F = 



dD 



h(n)dK- 



d\- 



oo 



}. 



(6-12) 

where all the objects are associated to curve Cq. Metric coefficient e 2k is, up to an 
arbitrary constant factor, given by the following expression: 



290 , , 

e 2k = {det Aq}~2 T\ \\j - £|-S9 ( & h( K )(U\J {K) - ejdn 
f=i V/aD 



B 



x exp 1 7riBn ( f h(K)dHc) +f f (W^a* \ K K * — 2 In \k — k\) dndk 1 , (6.13) 
I ^JdD ' JdDJdD " J 

where B is the matrix ofb-periods and dll is the normalized basis of holomorphic 1- forms 
on curve C 



Proof. Formulas ( 6.12| ) and ( |6.13| ) are direct continuous analogs of ( |6. 9| ) and ( |6.10| ), 
respectively. Term 21n|K — k\ is subtracted, using the freedom to renormalize e 2k with 
an arbitrary (^^-independent factor, to achieve convergence of the double integral at 
k = K. 
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